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ABSTRACT 

We have used Suprime-Cam on the Subaru Telescope to conduct a V- and /-band imaging survey of 
fields sampling the spheroid of the Andromeda galaxy along its south-east minor axis. Our photometric 
data are deep enough to resolve stars down to the red clump. Based on a large and reliable sample of 
red giant stars available from this deep wide-field imager, we have derived metallicity distributions vs. 
radius and a surface brightness profile over projected distances of R = 23-66 kpc from the galaxy's 
center. The metallicity distributions across this region shows a clear high mean metallicity and a 
broad distribution ([Fe/H]~ —0.6 ± 0.5), and indicates no metallicity gradient within our observed 
range. The surface brightness profile at R > 40 kpc is found to be flatter than previously thought. It 
is conceivable that this part of the halo samples as yet unidentified, metal-rich substructure. 
Subject headings: galaxies: individual (M31) — galaxies: halos — galaxies: structure 



1. INTRODUCTION 

Our understanding of how disk galaxies like our own 
or the Andromeda galaxy ( M31) formed plays a key role 
in n ear-field cosmology ([Freeman fc Bland-HawthorrJ 
l2002f ). because such nearby galaxies offer us detailed 
views of galactic structures through their resolved stars. 
In particular, the components of the extended low sur- 
face brightness stellar halo, old field stars and glob- 
ular clusters, provide invaluable information on early 
chemo-dynamical evolution of disk galaxies over the past 
> 10 Gyr (i.e., well before thin disk components ap- 
peared). The spatial distribution of galactic halos sug- 
gests that proto-galaxies were much larger than the sizes 
of currently bright galactic disks, possibly a result of 
hierarchical assembly of subgala ctic systems orbiting at 
larger radii (|Searle fc Zinnl [i978). The importance of ac- 
cretion and merging of small systems in the galaxy for- 
mation process is also evident in the spatial and kine- 
matic substructures in the Galaxy' s stellar halo (e.g., 
lYannv et al]l200l iHelmi et alJfl999h . 

The volume density profile of the Milky Way halo is 
characterized by an r -35 power-law of Galactocentric 
radius, based on dire ct counts of h alo tracers or their 
orbital motions fe.g.. lHarris I I1976T ). By contrast, An- 
dromeda's spheroid ap pears to have a complex structure . 
In their imaging study IPritchet fc van den Berghl Jl994) 
found that the surface brightness (SB) profile along the 
minor axis is characterized by a de Vaucouleurs i? 1 / 4 
law over projected distances R of 1" to 1?5 (few pc to 
20 kpc ) from the M31 center. Later, [Guhatha kurta et aT] 
(2005) reported, based on spectroscopic selection of red 
giant branch (RGB) candidates in several halo fields, that 
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the SB beyond R ~ 20 kpc shows a flatter profile with 
R~ 2 - 3 , and that t he halo may extend to R > 150 kpc. 
Ilrwin et al.l (|2005l ) also found, from their photometry of 
RGB stars over R = 20-55 kpc along the minor axis, a 
comparably flat SB profile that is fit by an R~ 2 3 power 
law or an exponential law with a scalelength of 14 kpc. 
It remains to be seen whether this flat portion of the halo 
profile is affected by substructure, and if it is dominated 
by either metal-rich or metal-poor populations. The lat- 
ter issue is especially important since spectroscopic stud- 
ies suggest that M31's inner spheroid/halo has a radial 
gradient in mean metallicity, whereby metal-poor stars 
dominate at R > 50 kpc (e.g.. iKalirai et al1 l2006). De- 
tailed studies based on large numbers of stars (say from 
wide-field imagers) are required to investigate the funda- 
mental nature of M31's halo along its minor axis. 

In this Letter, we report on our Subaru/Suprime-Cam 
wide-field photometric observations of M31's halo aimed 
at obtaining statistically significant SB and metallicity 
distribution (MD) profiles along the south-east (SE) mi- 
nor axis of the galaxy. Our imaging survey is optimized 
to extract halo profiles in a region that is well outside 
the bright disk component and beyond the inner spheroid 
with the i? 1 / 4 brightness law. 

2. OBSERVATIONS AND DATA REDUCTIONS 



Using the Suprime-Ca m ([M ivazaki et al. 2002|) imager 
on the Subaru Telescope (jive et a l. 2004), which consists 
of ten 2000 x 4000 CCDs with a resolution of CC'202 per 
pixel and covers a total field-of-view 34' x 27', we have 
carried out a wide- field imaging survey of M31's spheroid. 
Our targeted fields are located between 23 and 66 kpc 
from the M31 center (Fig. [J). 

During four nights in August 2004, we obtained images 
of seven spheroid fields (hereafter referred to as Fields 1- 
7) and a control field in Johnson V and Cousins I bands 
with typical seeing FWHM of 0'.'9 - 1'.'3. These observa- 
tions were obtained in non-photometric conditions. We 
carried out additional imaging of the same fields on a 
photometric night in August 2005 in order to calibrate 
the data. The data were reduced with the software 
package SDFRED, a useful pipeline d eveloped to opti- 
mally deal with Suprime-Cam images (lYagi et all 120021: 
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Fig. 1. — Locations of our Subaru/Suprime-Cam fields (rectan- 
gular areas), overlaid on the surface de nsity map of M3 1 's red RGB 
stars from the INT/WF C survey of Ferguson ct al. (2003) an d 
adapted from Figure 1 of llrwin et a.1.1 T2005T K Adjacent Suprimc- 
Cam fields overlap by about 25 % to ensure photometr ic calibra- 
tion. In addition, the fields used by Irwin ct al. ( 2005|) for their 
surface brightness profile are delineated by the dashed lines. 



lOuchi et al.l 120041 ). Pho tometry o f our co- added images 
was calibrated through lLandoltl (|1992l ) standards, cor- 
recting for an airmass term and a color term in each 
field. 

We then conducted PSF-fitting photometry u sing the 
IRAF version of the DAOPHOT-H. software ([Stetson! 
|1987[) . We adopted a 3 a detection threshold for the ini- 
tial object detection/photometry pass and iterated the 
PSF-fitting photometry twice with 5 a and 7 a detection 
thresholds, respectively, in order to account for blended 
stars. A stellar PSF template was derived from about 
100 bright, isolated stars per image. Finally, we merged 
two independent V- and /-band catalogs into a combined 
catalog using a 1-pixel matching radius. 

It is worth noting that the morphological segrega- 
tion method cannot distinguish between M31 halo stars 
and and other point sources such as Galactic foreground 
dwarf stars, and compact extragalactic objects. To sta- 
tistically remove these contaminations from the targeted 
fields, we adopted a control field located at a Galac- 
tic latitude comparable to that of our targeted fields, 
on the assumption that it has the same abundance of 
foreground and background objects as that in the M31 
spheroid fields. 

Reddening corrections are applied in e ach fi eld based 
on the extinction maps of ISchlegel et al.l (|1998f ). and the 



iDean et al.l (|1978l ) reddening law E(V - 1) = l.ME{B 
V) and Aj = 1.31E(V - I). We adopt the M31 Cepheid 
distance m odulus of (m - M)„ = 24.4 3 ± 0.06 (~ 770 
kpc) from iMadore &: Freedmanl (|1995l ). We also note 
that our analysis of the SB and MD profiles below is un- 
affected by photometric incompleteness: we target RGB 
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Fig. 2. — Histog r ams o f the metallicity distributions, in the 
ICarrett a & Gratton ( 1997) scale, for the SE minor axis fields of 



Field 5-841 stars 
[Fe/H]„ v ,,= -0.62±0.55 
[FVH] m „,= -0.5±0.4 



Field 7 - 902 stars 
[Fe/H] IV „=-0.51±0.44 
[Fe/H]„. a =-0.5±0.3 




Field 6 - 1036 stars 
[Fe/H].„.= -0.60±0.52 
lFVH]„„, = -0.5±0.3 




M31's spheroid. In the upper left corner of each panel, we show 
the name of the field, the number of stars used to derive the MD, 
the average metallicity ([Fe/H] avc ) together with the associated 
standard deviation, and the median metallicity together with the 
associated semiinterquartile interval. 

stars brighter than / ~ 22.5 (mag), while our data are 
more than 90% complete down to I ~ 23 (mag). As an 
example, even if we count stars down to the 50% com- 
pleteness for the V band, the hitherto missing faintest 
and reddest stars are included in the analysis but the SB 
and MD profiles are changed by less than a few percent. 

3. RESULTS 

3.1. Metallicity Distributions 

Taking advantage of the large stellar samples available 
with the wide field of view of Suprime-Cam, we derive 
MDs of RGB stars in M31's spheroid, by comparing them 
to RGB fiducials define d by Galactic g l obular clusters 
in the same manner as Bcllazz ini et al.l (|2003t ). For a 
reliable secure determination of MDs using these tem- 
plates, we select RGB stars with —3.8 < Mj < —2.0 and 
0.9 < (V — I)q < 4.0. These selection criteria allow us to 
remove a number of contaminants, such as M31 asymp- 
totic giant branch (AGB) stars, AGB bump stars, and 
young stars, along with foreground and background ob- 
jects. We perform an interpolation procedure to obtain 
metallicities of objects in the spheroid and control fields, 
and subtract the derived MD of the control field from 
that of the spheroid field in order to statistically remove 
the effects of the foreground and backgro und contami- 
nants . For more detailed description, see Tana ka et al.l 
(f2007h . 

Figure [2] shows the MDs in Fields 1 to 7 of M31's 
spheroid. The vertical error bars denote a nominal uncer- 
tainty in each metallicity bin, yielding from the Poisson 
errors equal to ± y iV(spheroid + control) . It is worth 
noting that because of a large number of the RGB stars 
available from our Suprime-Cam data, the MDs obtained 
here have significantly small errors and thus are expected 
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Fig. 3. — The average metallicity plotted as a function of distance 
from the M31 center. Filled red circles present the average values 
of metallicity in the seven targeted fields, while open circles denote 
those in each subficld inside a Suprime-Cam field (see text). 

to be most likely compared with the results of the pre- 
vious studies, whic h examined more inner parts of the 
spheroid than ours (IDurrell et alJfeOOli iBellazzini et al.l 
2003; IDurrell et al.ll2004l) . 

It follows that the MDs have a broad distribution rang- 
ing from [Fe/H] ~ — 3 to the near solar metallicity. The 
average metallicity ([Fe/H] avc ) with a standard deviation 
and median metallicity ([Fe/H] mc( j) with a quartile devi- 
ation are shown in each panel . In addition, our F ield 2 
includes the field observed bv IDurrell et al.l (|2004r ). and 
the MD of Field 2 is almost consistent with their MD. 

The most striking feature of the MDs shown in Fig. [2] 
is their overall similarity. For the sake of comparison 
of them, we show their average metallicities as a func- 
tion of distance from the M31 center in Figure [31 Seven 
filled red circles present the average values of metallicity 
in our targeted fields, while many open circles are those 
in the subfields, which are obtained by dividing a single 
Suprime-Cam field into six fields (4000 pixels x 3000 pix- 
els per divided field) . These subfields help us to examine 
the fine spatial variation in a Suprime-Cam field. Con- 
versely, a filled red circle reflects the mean value of metal- 
licities of six subfields. Vertical error bars show standard 
deviation of the mean. These plot suggests that the stel- 
lar content of the spheroid has no metallicity gradient 
within our observational range. This may be consistent 
with the idea of chaotic merging of sm all building blocks 
first proposed bv lSearle &; Zi nn (1978;) for understanding 
the lack of a metallicity gradient in the globular cluster 
system of the Milky Way halo. 

3.2. Surface Brightness Profiles 

We estimate the SB profiles of M31's spheroid along 
the SE minor axis using the RGB stars, which are 
brighter than the AGB bump (/ ~ 23.3 mag). After 
extracting the secure RGB stars on which are imposed 
the same magnitude and color criteria as in preceding 
section, we divide them into the following two groups 
using the RGB templ ates of Galactic globular clusters 
(Bcl lazzini et a l. 2003): Metal-Poor group (MP) defined 
as -2.16 <[Fe/H]< -1.11 and Metal-Rich group (MR) 
defined as -1.11 <[Fe/H]< 0.07. The values of -2.16 
and 0.07 are both ends of the templates, and the tem- 
plate having a value of — 1 . 1 1 is chosen because it nearly 
marks a discontinuity in the MDs obtained here. 

In contrast to the inner dense parts of the halo in M31, 
it is possible to resolve the individual stars in its outer 
parts, because of their sparse density and of the suffi- 
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Fig. 4. — The left panel shows the V- and /-band minor-axis 
profiles for RGB stars as derived from the MP and MR group, 
while the righ t panel shows thos e based on the selection method of 
RGB stars bv llrwin et al.l (I2005T ). The V-band and /-band profiles 
are illustrated with blue and red points, respectively, which are 
derived from star counts in the magnitude and color selection boxes 
as described in the text. The error bars reflect a combination of 
Poissonian and background uncertainties. The dashed lines show 
an exponential profile with a scalelength of s = 21.9 kpc for MP 
and s = 70.0 kpc for MR (left panel), and s = A6A kpc for blue 
RGB stars and s = 58.5 kpc for red RGB stars (right panel). The 
dot-dashed line in the right panel s hows an exponential profile with 
s = 13.7 kpc by Irwin ct al. (2005), for the sake of comparison. 

ciently close distance of M31 to us. Thus the SB profiles 
in spheroid's outer part are available by directly count- 
ing individual stars. The sources of noise in this method 
of estimating the SBs arise from the Poisson statistics 
from the finite number of stars observed and contami- 
nation from Galactic stars along the line of sight to the 
spheroid. The contribution from unrelated Galactic stars 
can be estimated by observing other nearby fields: we use 
our control field at ~ 90 kpc to this end. 

The left panel of Figure 2] shows the SB profiles of the 
MP (blue) and MR (red) group, which are used to derive 
the V- and /-band profiles, respectively. The dashed 
lines show an exponential SB profile, exp(— R/s), with 
a scale length s of 21.9 kpc for MP and 70.0 kpc for 
MR. We note that at large distances R, the number of 
RGB stars is just slightly larger than that of background 
objects, so we can determine only the upper limit for 
the most of their errors. It is remarkable that the SB in 
the MP group decreases more steeply with R than that 
in the MR group, suggesting that as above stated there 
may exist a metal-rich substructure at R = 50 ~ 60 kpc. 

The right panel of Fig. 2] sho ws the SB profiles b ased 
on the different star selection bv llrwin et al.l (|2005l ): us- 
ing I/-band and Gunn i-band systems, they selected blue 
RGB stars with 20.5 < i < 22.5 and 24.85 - 2.85(V" - 
i) < i < 26.85 - 2.85(1/ - i) and red RGB stars with 
21 < i < 22 and i > 26.85 - 2.85(V" - i), and derived 
the V- and /-band profiles from these categorized stars, 
respectively. We shift /-band magnitude criterion about 
0.02 mag toward the fainter side in order to match the 
Gunn i-band system. We note that our selection cri- 
teria for MP and MR are more advantageous than this 
fixed magnitude and color selection of blue and red RGB 
stars for investigating the effects of the MDs on the SBs. 
To compare with the blue RGB star count profile in 
llrwin et al.l (|2005f ). which is fitted to an exponential pro- 
file with s — 13.7 kpc, we shift it vertically to roughly 
overlap with our blue RGB surface brightness profile in 
the radial range 1?6 — 3?0. This result shows that stel- 
lar density at around 45 kpc and 60 kpc significantly ex- 
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ceeds a mean halo profile obtained bv llrwin et alj {2005). 
Hence, it is possible that two as-yet-unknown substruc- 
tures consist of bright metal-rich stars. 

We note that the details of the SB profiles obtained 
here are found to be pretty insensitive even if we adopt 
much fainter V- and /-band magnitude limits, say down 
to V = 25 mag and I = 24 mag. Thus, the current 
results are real and robust against the effects of any con- 
taminations over the fields; our deep Suprime-Cam data 
enable us to assess the minor effects of these errors. 

4. DISCUSSION AND CONCLUSIONS 

As fossil records of galaxy formation, M31's spheroid 
offers us a global perspective about its spatial structures 
as well as stellar populations. Based on a wide-field im- 
ager of Subaru, we have confirmed the pre vious results 
(|Guhathakurta et al.l2005Hlrwin et alll2005| ) that the SB 
profile in the SE minor axis is much flatter at R > 23 
kpc than the innermost part of the spheroid described 
by the i? 1 / 4 law. Our finding of no metallicity gradient 
in this outer part over R = 23 ~ 66 kpc implies that dis- 
sipationless processes, such as accretion and merging of 
collisionless stellar systems, may have been at work for 
the formation of this spheroidal part, thereby erasing any 
metallicity nonuniformity. We note here that our control 
field at R = 90 kpc of the SE mi nor axis, which shows no 
substructures (jlbata et al.ll2007l ). may still contain some 
spheroidal stars, because M31's spheroid ma y extend up 
to R ~ 150 kpc (|Guhathakurta et al.ll2005f) . However, 
even so, the relative MD and SB profiles obtained here 
remain unaltered. 

Intriguingly, our finding of the possible overdensity re- 
gions at around 45 and 60 kpc is in agreement with the 
recent work by llbata et ail (|2007f ) (having appeared in 
astro-ph while preparing for the submission of our pa- 
per), who also found the stream- like features at similar 
locations, referred to as stream D and C, respectively. 
Our metallicity analysis suggests that their stream D 
appears to be slightly more metal-poor than stream C. 
Also, the MDs of stream C (our Field 6 and 7) have 
a high- metallicity peak sim ilar to those of the Giant 
Stream (jTanaka et al. 2007). However, it is yet uncer- 
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tain if stream C is indeed related to the Giant Stream 
because of the lack of the information on the distance to 
stream C as well as its internal kinematics. Deeper pho- 
tometry down to horizontal-branch magnitude and/or 
multi-object spectroscopy of stars will be important for 
clarifying the origin of the newly detected substructures. 

Our results reported here are not to t ally i ncon- 
sistent with th e find ing by iKalirai et al.l (|2006l ) and 
IChapman et al.l (|2006[ L who argue for the presence of 
metallicity gradient obtained from their kinematical 
studies of the bright RGB stars using Keck/ DEIMOS. In- 
deed, only our Field 2 overlaps Field aO of IKalirai et alJ 
(2006), where both fields yield basically similar metal- 
licity, and other fields of Field 2 to 7 along the minor 
axis do not overlap their survey regions (see their Fig- 
ures 1 and 12). It is interesting to remark that their 
Field m6 at R ~ 87 kpc along the minor axis shows 
rather high metallicity of [Fe/H]~ —0.85, being similar 
to that of their Field aO as well as ours. In contrast, 
other fields deviated from the minor axis (such as their 
al3, al9, and bl5) show low metallicities, which basically 
give rise to the reported metallicity gradient. Thus, this 
implies that the MDs along the minor axis fields in con- 
cern may be somewhat different from other outer halo 
fields, i.e., more higher surface density and more metal- 
rich than previously thought as reported here. A tempt- 
ing view for these properties of our observed fields is that 
there exists a metal-rich substructure at R of 50-60 kpc, 
corresponding to a fa int tail of the giant stellar stream 
( Fer guson "eIaI1[2002T l or a part of another stream. More 
extensive studies by observing large halo areas of An- 
dromeda, e.g., using much wider HyperSuprime camera 
under construction, will be necessary to arrive at more 
decisive conclusions and thus obtain the accurate forma- 
tion picture of this typical disk galaxy. 
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